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Abstract 

We consider a random matrix model of interaction between a small n-level system, 
S, and its environment, a A-level heat reservoir, R. The interaction between S and R 
is modeled by a tensor product of a hxed nxn matrix and a N xN hermitian Gaussian 
random matrix. We show that under certain ”macroscopicity” conditions on R, the 
reduced density matrix of the system ps = is given by ~ exp{—/3Hs}, 

where Hs is the Hamiltonian of the isolated system. This holds for all strengths of 
the interaction and thus gives some justification for using pg' to describe some nano¬ 
systems, like biopolymers, in equilibrium with their environment |26] . Our results 
extend those obtained previously in [mEZ! for a special two-level system. 


1 Introduction 

The properties of a system S, in contact with a thermal reservoir R, is an old yet perennial 
problem of statistical mechanics. Writing the Hamiltonian of the composite system S' U i? as 

Hsur = Hs G) 1r + Is ® Hr + VsR (1.1) 

where Hr and Hr are the Hamiltonians of the system and of the reservoir and Vrr is the 
interaction between them, the canonical Gibbs density matrix of S' U i? is given by 

pfvjR = exp{-/3V5uR}/^sus- (1-2) 

The corresponding reduced density matrix of the system is 

Ps = T^rRpf^R = ex^{-0s}/Zs. (1.3) 

Here Hr is the ’’effective” Hamiltonian of S', which will have the form 

Hs = Hs + % (1.4) 

and Vr will in general depends on /S, Hr, Hr and Vrr unless Vrr is ’’negligible” and ps can 
be replaced by 

pf = exp{-/3Vs}/Zs. (1.5) 
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The use of fll.51) for the density matrix of S may be appropriate even when Vsr is not 
small, if the system S is macroscopic, e.g., a system in a large box A, and the interaction 
Vsr takes place only along the boundary dA. Then the calculations of the properties of 
the system far from the boundary are approximately independent of the interaction Vsr, 
becoming rigorously so when A —)■ cxd (there are exceptions when the system is at a hrst 
order phase transition). Our concern here is however with the case when S is small so the 
above considerations do not apply. 

Recent technological advances, making it possible to create and manipulate meso- and 
nanosystems, including biopolymers, colloidal particles, etc. have brought the problem of 
a micro-system in contact with a reservoir to the fore. Likewise, in recent studies of the 
foundation of quantum statistical mechanics 0 [3 El [la ini i23], where the approach to 
equilibrium is related to the system-reservoir entanglement the notion of an equilibrium 
state for non macroscopic systems also plays an important role. In all such cases the nature 
of the interaction Vsr between the system and the reservoir is clearly important and ps 
is not necessarily of the Gibbs form fll.Sp . A specihc example of this importance is the 
collapse transition in polymers, which depends strongly on the nature of solvent not just on 
its temperature [TO] . 

We note here that the distinction between the cases fll.Op and fll.Sp was made clearly 
by Jarzynski [T2] who showed that his equality between the work done on the system by 
changing a parameter in Hs from A to R is given by the difference of the free energies Fs 
given by log Zs in fll.Op evaluated at the values of the parameter A and B. This means that 
when Vsr is not negligible then this difference will depend on Hr and Vsr and not just on the 
temperature of the environment. This distinction is sometimes blurred in the literature 
both experimental and theoretical where various aspects of this problem are considered. 
Many of them go under the name of stochastic thermodynamics, where the equilibrium 
state of a nano-system in contact with an environment at temperature f3~^ is sometimes 
implicitly assumed to be described by the distribution fll.Sp . see, e.g. reviews [131 IMl 126] 
and references therein. In fact however these studies do not consider the whole microstate of 
the small system. Instead, one argues that the reservoir degrees of freedom as well as some 
’’internal” degrees of freedom of the small system change very rapidly compared to those 
observed, which we denote by X. They can therefore be assumed to be in thermal equilibrium 
and their effect on the evolution of X is then just described via some stochastic term (noise). 
This means in particular that the only effect of the reservoir on the equilibrium state of the 
X variables is to specify its temperature, similar to the situation for macroscopic systems. 
This seems reasonable in some cases (but not in situations like those in the example of the 
polymer described earlier). One may think then of the interaction term as fluctuating so 
rapidly that its detailed nature ’’washes out” as far as its effect on the equilibrium properties, 
or even the time evolution, of its slow variables are concerned. 

In this paper we consider a simple model of such a situation. We do this by making the 
interaction between the system and reservoir random. We obtain the Gibbs distribution for 
an arbitrary hnite-Ievel quantum system S using the frameworks of equilibrium statistical 
mechanics and assuming that the reservoir has certain macroscopic properties. 

Note that in this situation we can start, instead of fll.2p . which itself requires some 
justihcation [T21 [13], with the microcanonical distribution for the composite system at a 
hxed interval A around a macroscopic energy Erur, he., with 

Psur{Esur) = Xa{Hsur)/^sur{Esur), (1-6) 
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where xa is the indicator of the energy interval 


^ — {Esur — <5/2, Esur + ^/2), S « Esur- 


(1.7) 


and 


^sur{Esur) — Trsufl Xa{Hsur)- 


We then obtain ps of fll.dp with the inverse temperature given by [U] 


/5 


OSr 

OEr 


( 1 . 8 ) 


(1.9) 


where 

Sr = \ogVLR{ER) ( 1 - 10 ) 

is the entropy of the reservoir and Er is its energy (since for a reservoir which is very large 
compared to the size of system we have Sr ~ Ssvjr and Er ~ Es\jr,)- 

The rest of the paper is organized as follows. The model is described in Section 2, the 
results and discussion are given in Section 3. Sections 4 and 5 contain the proofs. The model 
is an extension of that introduced and studied in our works [min], where a simple case 
of a 2-level S was considered both in the equilibrium and non-equilibrium setting. Here we 
restrict ourselves to the equilibrium and leave dynamic considerations for a future work. 

According to the above, the proofs are given for the microcanonical distribution fll.bp 
- fll.lOp of the composite system. The results for the canonical distribution fll.2p of the 
composite prove to be a simple corollary of those for the microcanonical distribution and 
require a version of a standard (in statistical mechanics) saddle point argument. 

Note also that large random matrices have been widely used to model a variety of com¬ 
plex quantum systems, including heavy nuclei and atoms, mesoscopic particles, quantum 
networks, graphs etc., i.e., not necessarily macroscopic, systems HI El [22]. Thus, our set¬ 
ting could model a qubit or a microcluster in a mesoscopic (or even nanoscopic) particle, a 
quantum dot, small quantum network, etc. 


2 Model 


We describe now the model and quantities that will be considered. Let EIr be a iV x iV 
hermitian matrix, be its eigenvalues and 


N .c 

um{E) = N-^Y.^{E-E^), / 

i=i 


VRi{E)dE = 1 . 


( 2 . 1 ) 


be its density of states normalized to unity. We assume that un converges as ^ oo to a 
continuous density v in the sense that for any continuous and bounded function / we have: 


lim 

N—^OO 



f{E)uN{E)dE 



f{E)u{E)dE, 



v{E)dE = 


1 . 


( 2 . 2 ) 


Let Wr be a random hermitian matrix, distributed according to the Gaussian unitary 
invariant law given by probability density 

Zpexv{-T:rWll2}, (2.3) 


3 









where is the normalization constant. In other words, we assume that the entries of 
hhi? = ^kj = hhj'A: are complex and independent for 1 < j < A; < Gaussian 

random variables such that 

E{W,t} = E{W%} = 0 , = (1 ( 2 . 4 ) 

This is known as the Gaussian Unitary Ensemble (see e.g. [22], Section 1.1). 

Let also Hs and E 5 be arbitrary n x n hermitian matrices. We define the Hamiltonian 
of our composite system S' U i? as a random nN x nN matrix (cf. fll.ip l 

Hsur = Hs G) Iw + In G) Hji + S 5 0 Wr/N^^'^, (2.5) 

where 1 ; is the unit I x I matrix. 

It should be noted that our results are valid not only for the special Gaussian distributed 
interaction fl2.3p - fl2.4p . but for any real symmetric or hermitian random matrix in fl2.5p . 
whose entries Wjk, 1 < J < A; < iV are independent and satisfy fl2.4p . However, in this case 
the techniques are more involved requiring an extension of those of random matrix theory 
for so called Wigner matrices (see [22], Section 18.3). 

In addition, our results are also valid in the case, where Hr is random and independent of 
Wr for all N. In this case we have to assume that the sequence {Hr} of random matrices is 
defined for all iV —)■ cx) on the same probability space as the sequence {Wr}, is independent 
of {W r} and satisfies fl 2 . 2 p with probability 1 . 

Having Hr and Wr random, we can view the second term in fl2.5p as the Hamiltonian 
of a ’’typical” iV-level reservoir and the third term as a ’’typical” interaction between the 
system and its reservoir. It is worth mentioning that the notion of typicality has recently 
been used in the studies of the foundations of quantum statistical mechanics, including the 
form of reduced density matrix in equilibrium [7] E] [191 HD 122] . In addition, the randomness 
(frozen disorder) is a basic ingredient of the theory of disordered systems. Its successful 
and efficient use is justihed by establishing the selfaveraging property of the corresponding 
results, i.e., their validity for the overwhelming majority of realizations of randomness for 
macroscopically large systems, see e.g. [T 8 ] . 

In our case the selfaveraging property is valid in the limit N ^ 00 and is given by Result 
I below. This suggests that in our model the iV —)■ cx) limit plays the role of the macroscopic 
limits in statistical mechanics and condensed matter theory. Note however that in statistical 
mechanics the density of states of a macroscopic reservoir is multiplicative in its volume 
and/or in its number of degrees of freedom. This and the macroscopicity of the system lead 
to the Gibbs form fll.5p of its reduced density matrix, if the system-reservoir interaction is 
of short range and confined to the boundary of S [15]. 

If, however, the system is small, then, as is already noted, its reduced density matrix is 
not fll.5p in general even if the reservoir is macroscopically large. This is well understood in 
statistical mechanics and in principle in the stochastic thermodynamics community [niHS]. 
This is especially in the frameworks of dynamical approach where the problem was first 
studied by Bogolyubov [1] for a classical system consisting of a harmonic oscillator interacting 
linearly with the macroscopic reservoir of harmonic oscillators, and then in a number of 
interesting and rather general quantum models with macroscopic many-body reservoirs, see, 
e,g., [DID IS]- Taking into account that according to ( 12.1 p and ( 12 . 2 p the density of states 
of our reservoir is asymptotically additive (but not multiplicative) in N, we conclude that 
the Gibbs form of the reduced density matrix seems to be even less likely in our case than 
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in the spin-boson model and that one needs additional conditions and procedures to obtain 


the Gibbs distribution fll.Sp . We will discuss such conditions below. Here we only mention 


that an asymptotically additive in volume density of states is the case in the one-body 
approximation of solid state physics. 

According to (TB] the Gibbs distribution for the two-level (n = 2) version of our model 
can be obtained if we assume that the reservoir consists of large number 1 << J << 
of independent or weakly dependent parts, more precisely, that the normalized density of 
states z/^r of the reservoir has a special ” quasi-multiplicative” form fl3.17p . One can view this 
assumption as an effort to obtain the multiplicativity of the density of states of the reservoir 
on a scale which is intermediate (mesoscopic) between the microscopic and macroscopic 
scales. It is shown below that the same assumption on 1/^ allows one to obtain the Gibbs 
distribution fll.Sp for an arbitrary hnite dimensional S. 


3 Results 

We are interested in the asymptotic form as iV —)■ cx) of the reduced density matrix 


(AT) _ N ^Tr^ Xa{Hsur) 

~ iV-lTrXA(J?su«) 


(3.1) 


corresponding to the microcanonical distribution (II.6p - (II.8p of our model composite system 



Our hrst result is 

I. There exists a positive definite, trace one, non-random matrix ps such that for any 
fixed A of fll.711 we have with probability 1 



(3.2) 


The result is proved in Section 4. 

To describe the limiting density matrix ps, we start from the relation 


Xa{Hs[jr) 


(3,3) 


where Thsur resolution of identity (spectral projection) of the hermitian operator Hsur 

(12. 5 p corresponding to the spectral interval A of (II.7p . Denoting 


efHdX) = 


(3.4) 


we can write, using fll.7p . fl3.3p and fl3.ip . 



Tl's 4"'(A) 


(3.5) 


Thus, we have to hnd the limit 


es(A) = lim e^^^(A) 


(3.6) 


5 














of the n X n random hermitian matrix e^^(A) and then we can write (13.21) as 

^ e5(A) 

Trs es(A) 

To find 65 (A) of (13.61) we will use the resolvent 

G(z) := (Hsur - ^ 0 


(3.7) 


(3.8) 


of Hsur- Indeed, we have, by the spectral theorem for hermitian matrices. 


G{z) 


(^z\) 

A — 2 ; 


7 ^ 0 , 


(3.9) 


and we write here and below integrals without limits for those over the whole real axis. It 
follows then from (13. 4 p and (13.9p that 


gf\z) :=N-^TTnG{z) 


I 


es^\dX) 
X — z 


^z 7 ^ 0 . 


(3.10) 


Note that 65 is the matrix valued measure assuming values in n x n positive definite ma¬ 
trices and is uniformly bounded in N, since for any A C M the matrix being an 

orthogonal projection in nN space l-is ® ^Ri is of norm one and according to (13.4p 

|l4''>(A)||<n. 


Here ||z4|| is the standard matrix norm of a matrix A and we used the inequality 


Tr^HII <nA||H||, 


valid for any hermitian matrix A in 'H 5 (g) "H/j. 

Recall now that for any non-negative hnite measure m on the real axis one can define its 
Stieltjes transform 

7 ^ 0 . 

which is analytic for ‘Az 7 ^ 0 and such that 

Q's(.s) Q'z >0, Q'z 7 ^ 0. (3-12) 


The correspondence between non-negative measures and their Stieltjes transforms is one-to- 
one, in particular 

m(A) = lim — [ As{X + i5)dX. (3.13) 

< 5 — 70 + 71 J ^ 

Besides, the correspondence is continuous with respect to the weak convergence of measures 
(see e.g. (12.21) 1 and the uniform convergence of their Stieltjes transforms on a compact set 
of C \ R (see e.g. [22], Proposition 2.1.2). 

It is easy to extend (13.111) - (I3.13P to the matrix valued positive dehnite and bounded 
measures and their matrix valued Stieltjes transforms, whose examples are (13.4p and (I3.10p 
respectively. Thus, to prove (13.71) it suffices to prove that for a compact set of C\R the matrix 
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valued functions flS.lOp converge with probability 1 to a non-random limit on a compact set 
in C \ R. 

Correspondingly, we prove in Sections 4 our second result 

II. Let e^^^(A), Hs and S5 be defined by \3.tA) - Hd.im) . \2.1\) - \2.2\) 

and h2.^) . Then: 

(i) the limit 65 (A) of l\3.6\) exists with probability 1, hence 113.2\) and 
same probability; 

(ii) if fs is the matrix Stieltjes transform of es, i.e., 

e 5 (A) = lim — [ ‘2sfs{\ + iS)dX (3-14) 

5^0+ 71 Ja 

then fs is a unique solution of the matrix equation 

fs{z) = I {E + Hs-z- Tsfs{z)Ts)-^ n{E)dE (3.15) 

in the class of matrix valued functions analytic for 7 ^ 0 and such that (cf. h3.12\) ) 

{fs{z) - fs{z))/iz -z)>0, ^Zj^ 0, (3.16) 

and for any hermitian matrix A we write A > 0 if A is positive definite. 

Remark. The case n = 1 of the above assertion corresponds to a particular case of the 
deformed Gaussian Ensembles of random matrix theory (see Section 2, Theorem 2.2.1 

in particular). The case n = 2, Hs = sa^, s > 0, and S 5 = where ax and a^ are the 

corresponding Pauli matrices, was considered in [mini, while studying a random matrix 
model of quantum relaxation dynamics. 

The obtained limiting reduced density matrix ps (see fl3.2p and fl3.7p i is generally not of 
the Gibbs form fll.5p . We thus need additional assumptions on the structure of reservoir in 
order to obtain the Gibbs distribution in the frameworks of our model. In our previous work 
[16], where the case n = 2 was considered, it was assumed that the reservoir consists of a 
large number J ^ 00 of practically non-interacting ’’macroscopically inhnitesimal” but also 
sufficiently large parts (i.e., a kind of ’’coarse grained” structure of reservoir). This can be 
implemented by writing the density of states fl 2 . 2 p of reservoir as the convolution of J copies 
of a certain density q: 

u ■= pj = q*fi (3.17) 

The fact that we are going to consider the asymptotic regime J —)■ cxo after the limit N ^ 00 
can be interpreted as a formalization of the inequality determining our intermediate scale 

n « J « N, (3.18) 

or, denoting Uq = N/J the parameter characterizing the ’’size” of inhnitesimal parts, as the 
condition Ug » 1 . 

A simple example of the above is the Gaussian density 


{3.'?) hold with the 
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where £o is assumed to be of the order of magnitude of the characteristic energies of ” macro- 
scopically inhnitesimal” parts of R, and a is their energy spread. In this case the density of 
states (13.171) of the reservoir is also Gaussian 

The formula makes explicit one more property of our ’’macroscopic” reservoir: its charac¬ 
teristic energies are of order of Jeq with the spread while those of Hg and S 5 are 

independent of J, hence, much smaller. 

The Gaussian ”ansatz” (13.191) - (I3.20p was used in [161IH] • fhis paper we study the 
equilibrium properties of 1 12.51) with an arbitrary hnite n and nxn hermitian Hg and S 5 in 
(12.51) . As already noted, our results are valid for a wide and natural class of ” macroscopically 
inhnitesimal” part distributions q in (13.171) . see Section 5. In particular, one can mention 
the densities 

(ld = qi{£) = ,0<e<eo, (3.21) 

where qi is the density of states of the one dimensional harmonic chain and qd is the density of 
states of the d-dimensional harmonic cubic lattice. Different choices of q can model different 
types of reservoirs. For carrying out the proofs we shall require two conditions on q: 

(i) q decays superexponentially at — cxo; 

(ii) there exists 1 < a < 2 such that 

J q°‘{e)de < 00 . (3.22) 

The conditions seem fairly natural. Indeed, condition (i) requires a rather ’’thin” if any 
spectrum of negative energies of large absolute value, thereby it is closely related to the 
stability of the reservoir. Furthermore, since (I3.22p is valid by dehnition for a = 1 (recall 
that g is a probability density), condition (ii) requires a certain regularity of the density of 
states, also usually assumed in statistical mechanics. Note also that if q is zero for large 
negative energies, then one can use the polynomial decay of the Fourier transform of q as 
another condition of its regularity. 

Here are simple examples of the above. The hrst is the Gaussian density (I3.19p . for which 
the validity of the Gibbs distribution in the asymptotic regime (I3.18P for n = 2 was proved 
in [16]. Here a is any real number of (1,2]. The second example is the exponential density 
q{e) = \ where again a G (1, 2]. The third example is the density of states (I3.2ip of 

the simple cubic d-dimensional crystal. Here a G (1,2) for d = 1 and a G (1, 2] for d > 2. 

Viewing the density of states (13.171) as the partition function of the microcanonical en¬ 
semble of an ideal gas of J particles having each the partition function q one can introduce 
an analog of the entropy per particle 

s{e) := lim J~^ \oguj{Je). (3.23) 

J—>-oo 

The existence of the limit, its continuity and convexity can be proved by now standard 
argument of statistical mechanics (see e.g. [25] l20]). 

One can also introduce an analog of the inverse temperature (cf. (II.9p ') 


/3{e) = s'{e). 


( 3 . 24 ) 




























The corresponding quantities for the canonical ensemble of ’’macroscopically infinitesimal” 
parts of our reservoir are the analogs of the partition function and the free energy per particle 


Zj(/3):= j /(/3) -(/3J)-MogZj(/3) =-/3-MogV>(/3), (3.25) 

where 

= J e~^^q{e)de. (3.26) 

In writing the partition function fl3.25p we took into account that the Laplace transform of 
the J-fold convolution vj of fl3.17p is Jth power of the Laplace transform of q. Condition (i) 
above ensures that ijj is well dehned for all non-negative (3. 

We have also from fl3.23p and fl3.25p - fl3.26p . 

/3/(/3) = min{e/3 - s(e)} (3.27) 

£ 

and 

s(£) = mm he{l3), h^{(3) = efi - f{(3) (3.28) 

Note that f3f and s are convex. The convexity of (3f is immediate, since fl3.25p and fl3.26p 
imply 

iPf)” = - = -d, d:= J{e- eofQ{e)de > 0 (3.29) 

where ^ 

Q(£) = e-^^q{e) (^j e-^^'q{e')de'^ > 0 (3.30) 

and Eq is the first moment of Q. As for the convexity of s, it follows from the proof of fl3.23p 
[25l|2^. 

We have also from (13.251) - (I3.28p 

he(0) = 0, h;(0) =£-£, (3.31) 

where 

eq{e)de 

is the mean energy of the macroscopically inhnitesimal components of the reservoir. We will 
assume that 

£ > £, (3.32) 

since then the convexity of he and (13.3ip imply that the minimum s(£) of he in (I3.28p is 
positive, hence our temperature fl3.24p is positive as well 

/3(e) > 0, (3.33) 

although negative temperatures can also be considered in the frameworks of our model. 

It turns out, however, that the above standard formulas of statistical mechanics of macro¬ 
scopic systems are not accurate enough to obtain the Gibbs form fll.5p of the reduced density 
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matrix fl3.14p - fl3.15p of S in our model fl3.17p of reservoir. This is because of the ” logarith¬ 
mic accuracy” of the formulas, see e.g. fl3.23p . giving the large J leading term of \ogvj{Je), 
while we will need below the more accurate asymptotics formula: 

vj{Je) = /ij(£)(l + o(l)), J ^ oo, nj{e) := (3.34) 


This formula can be viewed as the Darwin-Fawler version of the equivalence of the micro- 
canonical and canonical ensembles in statistical mechanics (see e.g. m)- 

In fact, we do not need to know that ^j{e) is the asymptotic (13.341) of vj^Je) to prove 
our third result below, i.e., the validity of the Gibbs form of the energy distribution of the 
system in our model fl3.17l) of the reservoir. Instead, we will just use yUj of fl3.34p as an 
ansatz. Since, however, our proof is applicable to all hnite n > 1 and all n x n hermitian 
Hs and S 5 in fl2.5p and since the case n = 1, and Hs = S 5 = 0 corresponds to the reservoir 
itself, thus, in this case 'jj{E) = we obtain the asymptotic formula fl3.34p as the 

simplest case n = 1 of (13.361) . 

III. Under the conditions of validity of Results I and II above and for the density of 
states of reservoir given by ( [3. j7| ) with q satisfying conditions (i) and (ii), h5.I\) and H3.2^) 
in particular, we have: 

(i) the reduced density matrix of and ^3. 7| ) has a well defined limit as the width 6 
of the energy shell A of ([i. 7| )) tends to zero, i.e., A —)■ {E} and J is large enough: 


r 3j{E) 

Pj{E) := Inn ps = 7 ^- 

A^{E} Tis'yjiE) 


where 65 (A) is given by (13.6p : 
(ii) the limiting relation 


lim 7 j(Je)//ij(£) = e 

J—>-oo 


lj{.E) = lim 65 (A), 

A->{E} 


_ ^-pHs 


with pj given in \3.3Jif) , hence the Gibbs form H1.5\) of the limit (in view of H3.35\) ) 

lim pj{Je) = 

J—>-oo 


Trse-hHs ’ 

and (EH) - (EH- 


(3.35) 


(3.36) 


(3.37) 


in which fl := I3{e) is defined by (I3.23P 
The result is proved in Section 5. 

Remark. The above result concerns the microcanonical energy distribution (II. 6 p 


(3. 


of the composite system and the corresponding reduced density matrix (13.ip . The 
case of the canonical distribution (II.2p of the whole system and the corresponding reduced 
density matrix, i.e., the passage from (II.4p to (II.5p can be readily obtained from (I3.35P - 
(I3.37p . Indeed, we have from (II.2p . the spectral theorem for Hsur, (13.4p and (I3.35P that 
with probability 1 : 


lim 


Tr^e 


-h^SuR 


= lim 


je-^^es{dE) _ j e-^^^j{E)dE 


N^oo Tr 5 TrRe I^sur n^oo J e h^Trses(dE) f e h^TTs'yj(E)dE 


(3.38) 


Using then (I3.34p . (I3.36P and a simple saddle point argument for integrals on the r.h.s. of 
(I3.38P yield 


lim 




J^oo J e-h^TTs'yj(Js)p/(s)de 
i.e., the Gibbs distribution again. 


Trse-h^s ’ 


(3.39) 
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4 Selfaveraging Property and Limiting Reduced Den¬ 
sity Matrix 


In this section we prove Results 1 and II above. We note that the corresponding assertions 
as well as their proofs are generalizations of those for the deformed semicircle law (DSCL) 
of random matrix theory, see [22], Sections 2.2 and 18.3. Just as in the case of DSCL the 
passage from the positive dehnite matrix measures (positive measures in the case of DSCL) 
to their Stieltjes transforms (see (13. 5 p and (I3.10p ) reduces the proof of (13.ip - (13.2p to that 
of (I3.15P - (I3.16p . The latter facts will follow if we prove that the variance of all the entries 
of (I3.10p vanish fast enough as n —?■ oo and that the expectation of (I3.10p converges to a 
unique solution of (I3.15P as n —)■ oo. 

We will use the Greek indices varying from 1 to n to label the states of the systems and 
the Latin indices varying from 1 to iV to label the states of reservoir. 

Then we can write n x n and nN x nN matrices of (I3.10p and (13. 8 p as 






(4.1) 


and 


We will prove now the bound 


N 


Qa/si^) — y ^ Gaj,i3jiz)- 

1=1 


(4.2) 


Var{( 7 „^(z) := ^{\g^p{z)\^} - \E{g^p{z)}\^ < Cs{z)/N\ (4.3) 


and we write here and below Cs{z) for quantities which do not depend on N and are hnite 
for ^z 7 ^ 0 . 

To this end we view every gap as a function of the Gaussian random variables 
of (12.3p - (12.4p and use the Poincare inequality (see [22], Proposition 2.1.6), yielding 


N 

Var{^„^(z)} < ^ E 

j,k=l 



(4.4) 


The derivatives on the right can be found by using the resolvent identity for G 


dgapjz) 

dWjk 


< 


1 

Ar3/2 


n 

[Gsp{z)Ga'y{z))jk^^S 

7,5—1 



\(GsPz)G„-,(z)),,\^ Y, 


7,5—1 


7 ,<5-1 


1/2 


where Gap denotes the N x N matrix 

Gap{z) = {GajM^)}f,,=V 

The n X n matrix S 5 = {'^ap}^ p=i is dehned by (12.5p and the second line results from 
Schwarz’s inequality. Taking into account the bounds 

|| G (^)|| < \Gaj,pk{^)\ < ||G(^)|| < (4.5) 
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valid for the resolvent of any hermitian matrix, it is easy to find an analogous bound for the 
matrix Gap{z)\ 

||G„^(^)|| < (4.6) 

The above relations imply for the r.h.s. of fl4.4p 

1 ” 

Var{9„s(z)} < ^Tr^El 

7,(5—1 

and using now (14.61) and the inequality \TiA\ < A^| |^| | valid for any TV x TV matrix, we obtain 

tA 

i.e., the bound fl4.3p with Cs{z) = n^Tr 5 S|/|$> 2 ;|^. 

Denote 

/f (2) = /osW = EfesW}. (4.7) 

We will now prove that for any compact set K C C\R there is a subsequence 
of analytic matrix functions, which converges on K as —?■ cx) to a solution of fl3.15l) . 

Repeating almost literally the argument, which leads to Req. (2.2.8) of [22], we obtain 

E{G(^)} = [{Hs - ^sff\z)^s) ®1 r + 1s®Hr- zlsui?)”' + Rn{z), (4.8) 

where the nN x nN matrix Rn{z) admits the bound 

\\RNiz)\\<Csiz)/N\ (4.9) 

Applying to fl4.8p - fl4.9p the formula 


N 

Tr_R(p(A 5 ( 8 ) !/{ + I 5 ( 8 ) Br) = + 6 ;), 

(=1 


valid for any function y9, any n x n matrix Ar and any N x N hermitian matrix Bji with 
eigenvalues and using the dehnition fl2.ip of the density of states of the reservoir, we 

obtain for the expectation fl4.7l) 



i']si{E)dE 

E + HR-z-T.Rff\z)T.R 


+ rjR^z), 


where 

IM^)||<C5(^)/fv^|^>^|^ 

It follows from 04.71) and 04.Sp that 


( 4 . 10 ) 



Hence, for any compact set K C C\M there exists a subsequence of bounded 

analytic matrix functions, which converges uniformly on K to an analytic matrix function / 
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and passing to the limit N ^ oo in fld.lOp we obtain fl3.15l) for z & K. The validity of the 
equation for any G C\M follows from the analyticity of / in C\M. 

Let us show that fld.lSp is uniquely solvable in the class of analytic matrix functions 
fl3.16p . Choose 

/L C G C : > IlSsll^} (4.11) 

and assume that there are two different solutions f'g and f'g, i.e., max^gx > 

0. It follows then from fl3.15p and 03.161) that 

1 < W^sW" j dEuiEME -Hg-Z- EgfgEg)-^\\.\\{E-Hg-z- Egf'gEg)-% (4.12) 

In addition, 03.161) implies the bounds \\{E — Hg — z — S 5 / 5 S 5 )“^|| < and \\{E — 

Hg — z — (see 05.4p for details). This and 04.12p yield the inequality 

1 < which contradicts 04.111) . 

The unique solvability of 03.15P - 03.161) implies that the whole sequence {fg^^}N con¬ 
verges to the limit /, a unique solution of 03.15P - 03.161) everywhere in C\M. 

It remains to prove that this fact together with 04.31) imply the convergence with proba¬ 
bility 1 of the sequence {g^^^}N of random analytic matrix functions to the solution of 03.15P 

- (ICT . 

It follows from the Tchebyshev inequality and 04.31) that for any 5 > 0 and any non-real 

P{ll/C(2) -sfT)ll > ^} < Var{<,f>(z)}/i2 < Cs(z)/i2|3z|‘>V^ 

Hence, for any non-real z the series 

OO 

N=1 

converges for any <5 > 0, and by the Borel-Cantelli lemma and the convergence of 
to fs we have with probability 1 for any non-real 2 ; 

lim gf\z) = fs{z). (4.13) 

N—^oo 

Let us show that g^^^ converges to fs uniformly on any compact set of C\M with probability 
1. Because of the uniqueness of analytic continuation it suffices to prove that with the same 
probability the limiting relation limAr^oo g^^\^j) = fsi^j) is valid for all points of an inhnite 
sequence “^Zj > rjo > possessing a hnite accumulation point. Denote r2(.2) the 

set of realizations, where fl4.13p is valid. According to the above P{r2(zj)} = 1, Vj. Hence, 

j>i 

This proves the uniform convergence of g^^"' to fs on any compact set of C\M with probability 

1 . 

It is known that the one-to-one correspondence between the non negative measures and 
the Stieltjes transforms is continuous in the topology of vague convergence of measures and 
the uniform convergence on a compact set of C \ M [22], Proposition 2.1.2. An analogous 
assertion can be easily proved for the positive dehnite matrix measures and their Stieltjes 
transforms. By using this fact and the uniform convergence with probability 1 of g^^'^ to fs 
proved above, we obtain 03.61) . 
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5 Gibbs Distribution 


We will prove here our Result III, i.e., relation fl3.37p . by using Result II (see fl3.15p - (13.1 bp ) 
and our model 03.171) of the reservoir. Recall that we assume that the density q of 03.17P 
satishes the conditions (i) and (ii) of Section 3, (15. 2 p - (I3.22p in particular. 

We will now use the matrix identity (see formulas (15.4p for its justihcation) 

— = i > 0, 

C Jo 

with C, = z + Tjsfs{.z)Tjs — E — Hs to write the basic equation 03.151) with > 0 as 

poo _ ^ 

fs{z) = / H = Hs- ^sfs{z)Es, (5.1) 

Jo 

where 

= j e~^^^q{e)de, (p(0) = 1, |9?(t)| < 1, 

and we took into account that the Fourier transform of the convolution 03.171) is 
proofs in this section are based on the above integral representation of fs- 
Let us prove hrst that the limit 

fs{E):= lim fs{E + z6) (5.3) 

<5->-0+ 

exists and is bounded for all real 77 if J is large enough. To this end we hrst use the bounds 

||(M - C}-'|| < (3C)"\ SC > 0; ||e-**"^|| < 1, t > 0, (5.4) 

valid for any complex matrix such that ‘JsM > 0. Indeed, in this case > 0. Matrices 

with this property are known as accretive and for them the corresponding bounds \ \{—iM + 
C}~^ll < (3f^C)~^) 3^C > 0; ||e“*^|| < 1, t > 0 are valid (see e.g. [H], Section IX.1.6). Noting 
that ‘Jsfs{z) >0, Q'z > 0 implies ‘Js'Efs{z)T, > 0 for > 0 and using the second bound of 
05.41) . we obtain from 05.11) 


(5.2) 

Our 


ll/s(j)ll < / Mt'pdt, Sz > 0. 


(5.5) 


It follows then from 03.22p and the Hausdorff-Young inequality that 

\ Jo/a- 


( f* \ '70/O' 

/ q°‘{e)de] < oo, Jo = a/{a — l) 


and we obtain from 05.2p . 05.5p and 05.61) 


(5.6) 


1 1 / 5 ( 2 ;)11 < / \(p{t)\-^°dt < 00 , J > Jq. 

Jo 


(5.7) 


It is important that the r.h.s. of the above bound is independent of .2 for > 0 and J > Jo, 
thus the same bound holds for fs{E) of 05.3p with any real E and uniformly in J > Jq. 
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We will use now a stronger version of the inversion formula fl3.13p corresponding to 
measures with a bounded density m': 

m'(A) := lim m(A)/|A| = lim 7 r“^ 9 's(A + Z(5). (5.8) 

Since the same formula holds also in the matrix case, we will divide the numerator and 
denominator of fl3.7p by |A| and pass to the limit A —)■ {E} to obtain fl3.35p . 

In fact, an elaboration of the above argument allows us to proves the relation 

lim ||/s(Je)|| = 0, (5.9) 

J—>-oo 

valid for the limit fl5.3p with any e varying over a hnite interval which will be used below. 
Indeed, it is shown in the Appendix that if satishes fl3.22p . then for any to > 0 

max |v5(t)| = (po < 1- (5.10) 

t>to>0 

Choose a 5 > 0 and write the integral in fl5.5p as the sum of integrals over (0, 5) and {6, oo). 
Then it follows from fl5.2p . I3.22p and fl5.6p that the hrst integral is bounded by 6 and the 
second is bounded by 

noo 

^0''° J > Jo, (5.11) 

where C is independent of J and 8. Thus, passing to the limit J ^ oo and then 5 —)■ 0, we 
obtain fl5.9p . 

In fact, a bit more careful calculation yields that the r.h.s. of fl5.9p is of the order J“^/^ 
if the hrst moment of q is zero, and J~^ if the hrst moment is not zero. This can be seen 
also in examples (13.19p and fl3.2ip . 

Introduce the real and imaginary part of fs{J^) (note that fs{J^) is well dehned in view 
of dES])); 

fs{Je) = n + tI (5.12) 

and take into account that according to our assumptions the function ip of fl5.2p can be 
analytically continued into the lower half-plane in t. This allows us to write the integral for 
X, which determines yj of fl3.35p according to (15.8p . as the sum 

I:=^fs{Je) = Ii + l 2 (5.13) 

of integrals over ( 0 , —i(3{e)) and (—ir/3(e), —i(3{e) + cxo), where /3(e) > 0 is the point where 
the function of (I3.28p achieves its minimum, see (I3.3ip - (I3.33p : 

Changing t to —ir in Ji we obtain 

rti 

I^= (5.14) 

Jo 

where we write 

e~'^^ = C{t) + iS{T). (5.15) 

for H of (15.ip and (3 instead of /3(e). Write also 

H = R -\- il, R = Hs — S 572 .S 5 , I = S 5 XS 5 > 0. (5.16) 
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If H were a complex number but not a matrix (e.g., if n = 1), then 


S{t) = - e --^ 

and 


h = - 

This and the inequality (see (13.281) - 


sin rJ = —e c{t), c(r) < r 

(5.17) 

rP 

■ / c(r)dr. 

Jo 

(5.18) 

(1332])) 


h(r) > (e - e)T, 

(5.19) 


would imply the bound 



(5.20) 


Comparing this with fl5.13p . we would conclude that the contribution of Ji in X of fl5.13p is 
negligible as J —)■ cxo, i.e., 

X = I 2 (1 + J —^ (5-21) 

We are going now to hnd a matrix analog of the above asymptotic relation. 

We note hrst that according to Appendix the matrix analog of fl5.17p is 

5(r) = - ||C(s)|| < se"ll^llcosh(s||/||), (5.22) 

Jo 


hence the matrix analog of fl5.18p is in view of fl5.16p 


h = - f e-^^-^^^IC{s)ds 

Jo Jo 

= f [ e-^^-^^^XsXJ:sC{s)ds. 

Jo Jo 

It is convenient to view the r.h.s. of the above formula as the result of application to the 
nxn matrix X of the linear operator A acting in the space 7V1„(C) of complex nxn matrices, 
i.e., to write the formula as 

h = A(X). (5.23) 

Then we have from (I5.22p just as in the scalar case (I5.17p - (I5.20p : 

llJill = ||A(X)|| <e^ll^ll||S5in|X||cosh(/3||S5in|X||)/((£-£)J)2. 


By using the bounds ||i7|| < 2||i7|| and ||X|| < 1, which following from fl5.9p . we conclude 
that the matrix analog of bound fl5.2ip is also valid, i.e., the contribution of A to the r.h.s. 
of fl5.13l) is negligible as J —)■ cxo in the matrix case as well. This fact can be expressed via 
the operator A: 

||A||<C'5/((£-£)J)^. (5.24) 

Consider now I 2 of fl5.13p . Changing the variable to t = —itq + a and using /ij of fl3.34p . we 
obtain 

_ _ poo _ 

l2/^i.j{e) = ^/2'Ks"{e)me-^^ / (5.25) 

Jo 
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where 


X{a) = isT + log'iljfj^a), '0^(a) = '0(/3 + *a)/'0(/3). (5.26) 

We write the integral on the right of fl5.25p as the sum of integrals over (0, ao) and (ao, oo), 
where ao is small enough; 

^2 — I 2 I + -^22- (5.27) 

In the hrst integral we use the expansion x(<^) = + o(l)), a —)■ 0 of (I5.26h 

yielding 

I 21 = 7r3?e“^'^(l + 0 ( 1 )), J —)■ cxo, 
and then fl5.9p and fl5.16p imply 

hi = 7re“^’^®(l + 0 ( 1 )), J —)■ 00 . (5.28) 


To deal with the second term 


1,2 = ^ 


'(TO 


of fl5.27p we use first fl5.9p and fl5.16p to obtain the bound ||e < K{e), where K{e) 

is independent of J. This and fl5.25p - fl5.26l) yield 


|^ 22 //^ j (£)|| < ^/ 271 shhJKh) / da, 


(5.29) 


where 


hfiih = hifd + ih/h{fd)- 

To estimate the integral on the right of fl5.29p we will follow the scheme of proof of fl5.1ip . 
which was based on the bounds fl5.10p and fl3.22p . Thus, we have to prove the analogs of 
these bounds for ifip. To this end we use fl3.26p and fl3.3Up to write 


h/sih ^ J ^ ^^'^Q{e)de. (5.30) 

We conclude that the ipjs is the Fourier transform of a non-negative function of unit integral, 
thus it satishes fl5.2p . i.e., |'0/3(o')| < 1, "^^(0) = 1. It follows then from Appendix that in 
this case we have an analog of (15.101) as well: 

max \'iljp{a)\ := Ki < 1. (5.31) 

cr>(To>0 

Furthermore, according to fl5.2p and fl5.30p . an analog of fl3.22p for '0^ is 


jQ^(e)de = je-‘^‘,‘(e)de<oo. 

It is easy to see this is indeed true in view of conditions (i) and (ii) of Section 3, (I3.22p in 
particular. 

This, fl5.27p and fl5.28p yield h = 7re“^-^s(l -|- o(l)), J —)■ cx). Now use fl5.23p to write 
fl5.13p as X =A(X) - I - J 2 . Then fl5.24p and fl5.25p imply 

X//ij(£) = (1 - Ayh 2 /fij{e) = 7re"^^®(l -h o(l)), J ^ 00 . 

This and fl5.8p yield fl3.36p . It remains then to divide the numerator and the denominator 
in fl3.35p by fij and to use fl3.36p to obtain (I3.37p . 
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Appendix 


(i) Proof of We prove here that if q is non-negative, continuous and of unit integral, 

then we have flS.lOp in addition to ( I5.2p . 

Indeed, the equality = 0 for some f ^ 0 implies the equality 

J~ cos et')q{e)de = 0, 


which is impossible for any non-negative and non-zero q satisfying fl3.22p . This implies fIS.lOp . 


(ii) Proof of 15. 2^) . The Duhamel formula 




valid for any two (generally non-commuting) matrices, yields the formulas for the terms of 
the r.h.s. of flS.lSp 

5(r) = - [ e-^^-^^^IC{s)ds, 


and 


hence 


C{t) = + / e-(^-*)^J5(s)ds 


C(s) = ' 

Jo Jo 


'•^ 21-1 




These formulas combined with the standard upper bounds for the terms of the series lead 
to (Km . 

In the commutative case the formulas are 


e sinrJ = e ^^/c(r), c(r) = / cossJ ds, 

Jo 

i.e., coincide with (I5.17p ). The non-commutative analog (I5.22p of (I5.17p is more rough, since 
we do not take into account the alternating signs in the above series, hence the effects of 
strong cancelations of its terms. 
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